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The upper critical field H C 2(T) of the tight-binding electrons in the three-dimensional lattice 
is investigated. The electrons make Cooper pairs between the eigenstates with the same energy 
in the strong magnetic field. The transition lines in the quasi-one dimensional case are shown 
to deviate from the previously obtained results where the hopping matrix elements along the 
magnetic field are neglected. In the absence of the Pauli pair breaking the transition temperature 
T C (H) of the quasi- two dimensional electrons is obtained to oscillationally increase as the mag- 
netic field becomes large and reaches to T c (0) in the strong field as in the quasi-one dimensional 
case. 
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In the semiclassical approximation, the upper criti- 
cal field H C 2(T) of isotropic superconductors is derived 
from Ginzburg-Landau-Abrikosoy-Gor'kov (GLAG) the- 
ory.ErB Lawrence and DoniachEP have applied GLAG 
theory to layered superconductors. Klcmm et al.u and 
Bulaevskii and GuseinovEP have shown that -ff c2 (T) dis- 
plays upward curvature when the magnetic field is ap- 
plied parallel to the layer. When the coherence length 
£(T) of the order parameter becomes smaller than the 
distance of layers, H C 2{T) is infinite. However, these 
results are based on a semiclassical approximation, in 
which only the phase of the wave function of electrons is 
changed by magnetic field. 

It was shown that the quantum effect of the electrons 
in the BCS thepxv. lead, tp reentrant behavior in high 
magnetic field. BQ&BEIllIJ The reentrance of the super- 
conductivity due to-|Landau quantization has attracted 
theoretical interest BB In that case, Cooper pairs are 
formed between electrons at the lowest Landau level in 
strong magnetic field. However, when only the lowest 
Landau level is filled in the three dimensional case, the 
system can be treated as a ID system, i.e., energy de- 
pends only on the momentum parallel to the magnetic 
field, and it has been shown that the system is unstable 
to thq-dposity wave state rather than the superconduc- 
tivityMe 

Recently, superconductivity in strong magnetic field 

Q observed in organic superconductors (TMTSF) 2 PF 6 . 
'cB Similar r^fuft has also been observed in 
(TMTSF) 2 C10 4 .BEJ Organic superconductors (TMTSF) 
are well described by quarter-filled tight-binding elec- 
trons (actually holes) with t a ~ 3000K, tb/t a ~ 0.1, 
tc/t a ~ 0.003. Since the magnetic field of 10T corre- 
sponds to 0/0o ~ 1/1000 in (TMTSF) 2 X, where 4> is 
flux per unit area and 4>q is the flux quantum, the ef- 
fect of the Landau level quantization is negligible. How- 
ever, when the magnetic field is applied perpendicular to 
plane, field-induced spin density wave (FISDW) is sta- 



bilized. The existence of FISDW shows that the quan- 
tum effect is important in quasi-one dimensional (Q1D) 
systems and the semiclassical approximation of the mag- 
netic field is not appropriate in these systems. 

LebecP has predicted that the Q1D superconductors 
should exhibit superconductivity in a strong magnetic 
field. The superconductivity observed in (TMTSF) 2 PF 6 
and (TMTSF) 2 C10 4 is thought to be the realization of 
the Q1D superconductivity in strong magnetic field. The 
mechanism of the reentrance of Q1D superconductivity 
is similar to that of FISDW. In the presence of the mag- 
netic field, the dimensionality of the system is reduced. 
When the magnetic field is applied in the c direction in 
the system with hopping matrix elements t a 3> tb 3> t c , 
the effect of tb, which makes the nesting of the Fermi sur- 
face imperfect with the imperfectness parameter of the 
order of t\/t a , disappears and spin-density- wave is in- 
duced. When the magnetic field is applied in b direction, 
the nesting of the Fermi surface stays imperfect but the 
orbital frustration is removed if we take account of the 
eigenstates in the magnetic field. In the Q1D case, the 
magnetic field necessary for the reentrance of the super- 
conductivity is much smaller than thaLin the case of Lan- 
dau level quantization. Dupuis et aZ.EP have extensively 
studied the mean- field transition line T C (H) of Q1D su- 
perconductors and shown the cascade transitions to the 
superconducting states. We have studied the anisotropic 
superconductivity with line nodes of the energy gap, i$ 
which is thought to be realized in orgaiiic superconduc- 
Xors (TMTSF) 2 X. In these papers,El&B) however, only 
warping of the Fermi surface in k z direction normal to 
H has been considered, and that in k y direction is ne- 
glected. 

As tb increases, quasi-one dimensional system turns 
into quasi-two dimensional system. Quasi-two- 
dimensional (Q2D) electrons, for example, /3-(BEDT- 
TTF) 2 I 3 (t a ~ 3000K, t b /ta ~ 0.5 and t c /t a ~ 0.02) 
have a two-dimensional cylindrical Fermi surface with 
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a weak warping along k z direction. When the magnetic 
field H is applied in b direction, the electrons move either 
in open orbit reaching the zone boundary or in closed or- 
bit as is schematically shown in Fig.l. 




Fig. 1. Fermi surface of a quasi- two dimensional conductor is 
drawn. In the presence of magnetic field along k y , the electrons 
move open orbit (a) and closed orbit (b). 



ta >tb ^ t c . We neglect the held dependence of /x, since 
the energy gap due to the magnetic held is very small 
in the anisotropic system with 4>/4>o *C 1 except for the 
bottom and the top of the band. 

In this paper the magnetic held H is applied in the 
b direction. We take the vector potential A as A = 
(0, 0, —Hx) and the non-interacting Hamiltonian is writ- 
ten as 



Ho 



E^i 

a.k 



V* 



M_i V 
V* Mq V 
V* Mi 



V v 



Ca, (3) 



where 



M n = —2t a cos[a(k x + nG)] — 24 cos(bk y ) — fj, — afisH, 

(4) 



V = -t r e h 



ct = (..., c t(fe-G),ct(fe) )C t(fe + G) 



(5) 



(6) 



The former gives the similar effect as in Q1D supercon- 
ductors. It have been discussed by LebectP and Dupuis 
et al. a that the Q2D superconductors will evolve from 
the GLAG .region to reentrant phase. Recently, Lebed 
and YamajiliiP have calculated the mean held transition 
temperature of Q2D superconductor and shown the reen- 
trant behavior. They have used the parabolic band in the 
conducting plane, i.e. the lattice structure is neglected 
in the plane. 

In this paper, we calculate the mean held transi- 
tion temperature numerically by taking account of the 
eigenstates of three-dimensional tight-binding electrons 
in magnetic held. In this formation, we can treat both 
Q1D and Q2D systems by changing tb/t a . 

We start from the tight-binding Hamiltonian (we take 
h = fes = Co = 1 , where Co is a light velocity, throughout 
the paper), 



-t c 

(i,j)o,c 



i.er 



(1) 



where c\ a and ci, a are creation and annihilation opera- 
tors, [L is chemical potential and cx/iB-ff is the Zeeman 
energy for \ (j) spin (a = +(— )) and 



4>o 



(2) 



In the above A is the vector potential. We consider 
the anisotropic system with the hopping matrix elements 



G = (G,0,0) 



2tt 

a 0o 



-,0,0 



(7) 



The creation operators of electrons can be written 
in terms of the creation operators of the eigenstates 
(*I(n,fc)) of eq.(3) as 

C t(fe + mG) = e imck * £ <j>* kx>K (m, n)*t(», (8) 

n 

where m and n are integers. 

Using eq.(8), the real space one-particle Green's func- 
tion is given by 

&.,*.(m,n)0t k ,(m',n) 



= EE 

k,n m > m ' 



1U} 1 - £ n.k.a 



x e 



i(m' -m)ck z i(r' ~r)-k+i(m'r' -mr)-G 



(9) 



where u>i = (21 + l)itT is a Matsubara frequency and I is 
integer. In this paper, we neglect the Zeeman energy for 
simplicity. The coefficients 4>k as .k z { , m,n) and eigenvalues 

£ n ,k = e ( n > fc x> k z ) - 2t b cos(bk y ) - fi, (10) 

can be calculated by diagonalizing the matrix in eq.(3) 
numerically, where e(n, k x ,k z ) is the eigenvalue of eq.(3) 
for — 2t(, cos(bk y ) — fj, = 0. If (j>/(f>o = p/q, where p and q 
are mutually prime integers, the matrix size of eq.(3) 
is q x q and the magnetic Brillouin zone is given by 
-n/(qa) < k x < n/(qa). 

Since we are interested in the instability to the super- 
conductivity in the quarter-filled tight-binding electrons 
at the temperature much smaller than the band width, 
the eigenstates near the top of the band are not impor- 
tant in the anisotropic hopping case (t c <C t a ). Thus, we 



Superconductivity of Ql and Q2 Dimensional Tight-Binding Electrons in Magnetic Field 



3 



calculate only 3/4 ~ 1/2 of the states from the bottom of 
the band in eq.(3) by neglecting V and V* for c(k + nG) 
with \k x + nG\ > (3/4)?r/a or \k x + nG\ > (l/2)n/a. 
We have checked that the result is not changed by this 
approximation. In this approximation, we can take the 
matrix size as ~ [tt/G] x [tt/G] and — G/2 < k x < G/2. 
The eigenvalues and coefficients do not depend on k z in 
this approximation and we write e(n,k x ,k z ) = e(n,k x ) 
and 4>k x ,k z (to, n) — </>fc x (m,n) (4>k x (m, n) can be taken 
real). 

In the mean field approximation, the linearized gap 
equation for s-wave pairing in coordinate representation 
is obtained as 

A(r) = X J dr'K(r,r')A(r r ), 

where 

K(r, r') = Tj^Gir, r', m)G(r, r', -iw,), 



(11) 



and A is the coupling constant. We write 



A(x) 



e^ +NG ^A N (q x ), 



(13) 



where q x is taken as —G/2 < q x < G/2. For even N, 
using eqs.(9) and (13), linearized gap equation is written 
as a matrix equation 



j' 



(14) 



where 



■ j,n)4> kx (m- j',n) 



k x ,k y n,n' rn 

x <\>-k x (~ m - 3, n ')<t>-k x (-m - j', ra') 

1 — f \^-n,k x ,k y ) ~~ f\^n', — k x — q x , — ky) 



X 



2(£n,k x ,k y ~\~ E-n' , — k x — q x 

where /(e) is the Fermi distribution function. 
For odd N, we get 

^■2j+i{q-x) = U 2 j + i } 2j'+iA 2 j'+i(q x ), 



(15) 



(16) 



where 



we can take the states (n,k) and (n', — k) which have 
the same energy e n ,k x ,k y = £n>,-k x ,-k y - Therefore the 
superconductivity is not destroyed by the orbital frus- 
tration in the strong magnetic field. This is the similar 
mechanism for the FISDW.tIP In the weak magnetic field 
the coefficient <f>k x (m, n) becomes small and the present 
results reproduce the GLAG results. 

For the Q1D superconductors with open Fermi surface, 
the energy dispersion in k x is taken to be linear and 
the Fermi velocity and the density of stateSpare taken to 
be constant in the previous calculations.&LrE£f As a first 
approximation, we take 



sgn(n)v F (k y )(\k x + nG\ 



(18) 



where VF{k y ) = 2t a a sin akp(k y ) and kp(k y ) are the 
Fermi velocity and Fermi wave number depending on 
k y , respectively. We may diagonalize the matrix for 
(12) k x + nG ~ kp(k y ) and k x +nG ~ —k-p(k y ) independently, 
as in the previous calculations. Then the eigenstates are 
given with the coefficients, 



n 



2j- 



k xy k y n,n' m 

0-fc*(- m - 3 - ±,n')<P-k x (-m - j' - l,n') 

1 — f( £ ii,k x ,ky) ~ f{£n', — k x —q x ,-k v ) 
2{Sn,k x ,k v ~t~ ^ri ,~k x — q x . — ky) 



The transition line is given by 1 — g\ = for eq.(14) 
and eq.(16), where g is the maximum eigenvalue of the 
matrix II of the even part or the odd part. In this paper, 
wc calculate the field dependence of the effective coupling 
constant g at low temperature instead of calculating the 
transition temperature. In the BCS theory Cooper pairs 
are formed by electrons with wave numbers k and — k 
and the energy of these states are different in the pres- 
ence of magnetic field, which causes the orbital frustra- 
tion. In the formulation of eqs.(14) and (16), however, 



4> kx {m,n) = J„ 



2t P . 



v F (k y )G 



and 



4> kx (m, n) = J- 



m-\-n 



2t c 



VF(k y )G 



for n, m > (19) 



for n, m < 0, (20) 



where J is the Bessel function. Note that in this approx- 
imation cf>k x (m, n) does not depend on k x but it depends 
on k y through vp(k y ). In Fig. 2, we plot the effective cou- 
pling constant g/go obtained by this linearized dispersion 
by solid, dashed, and dot-dashed lines as a function of 
aG/(27r), where go is the effective coupling constant for 
t c = , which corresponds to that in the absence of mag- 
netic field. 



— t b / t a = 0.0 (approx) 
-- t B /t, = 0.1 (approx) 

- t„ / t a = 0.3 (approx) 
O t b / 1, = 0.0 
41(1=0.1 



-i,2i'+i(?x) = fa* ( m ~ n ^fa* ( m ~ j'' 




(17) 



Fig. 2. Effective coupling constant as a function of aG/(2n) in 
the case of t c /t a = 0.02 and T/t a = 0.001. In strong magnetic 
field, the effective coupling constant obtained by diagonalizing 
the even part of the matrix II reaches to that of zero magnetic 
field. The coupling constant for odd part gives zero in strong 
magnetic field. Solid, dashed and dot-dashed lines are obtained 
in the approximation of linearized energy dispersion. Circle, tri- 
angle and diamond symbols are obtained without that approxi- 
mation. 



In strong magnetic field, the effective coupling constant 
of the even part is increased and that of the odd part is 
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decreased for each tp/t a - In the case of 4 = 0, we get 
the previous result.Ef As can be seen in Fig. 2, for larger 
4 / 4 the oscillation becomes small. 

Next, we calculate the effective coupling constant g/g 
by numerically diagonalizing the lower 3/4 or 1/2 of 
the matrix in eq.(3) without using the approximation 
eq.(18) and we plot the results as circles, triangles and 
diamonds in Fig. 2. For 4/4 = 0.1 the results are al- 
most same as that obtained by the approximation with 
the fc y -dependent Fermi velocity (eq.(18)) as expected, 
but the deviation is large for larger 4/4- 

We also study the quasi-two dimensional superconduc- 
tor with 4/4 = 0.5 and 1.0. In Fig. 3 we plot the effective 
coupling constant as a function of aG/(2ir). As is seen 
in Fig. 3, the effective coupling constant reaches to that 
for 4 = as magnetic field is increased. 



□ □ L / L = 1 .0 




simplicity. However, we can calculate the transition line 
with Zeeman term in this expression. The Zeeman term 
does not play any important role for the equal-spin- 
pairing case of the spin triplet. If the Zeeman energy 
is taken into account, the transition temperature of spin 
singlet is reduced due to the effect of Pauli pair-breaking 
except for the Q1D case (4 = 0). The superconductivity 
of Q1D systems is not completely destroyed, BBEiP since 
half of the density of states is available to make Cooper 
pairs _fbr_jthe Larkin-Ovchinnikov-F^ldfiJjerrell (LOFF) 
stateH^ as in the pure ID case-BBB 

In conclusion, we have shown the transition lines of 
quasi-one dimensional and quasi-two dimensional super- 
conductors in tight-binding model. The Green function 
in Q1D systems is described by the Bessel function if 
we apply the approximation that the energy dispersion 
in k x direction is taken to be linear. In this paper, we 
have shown that the Green function of the tight-binding 
electrons is numerically calculated without using the lin- 
earization of the energy dispersion. In the strong mag- 
netic field Cooper pairs are formed in the eigenstates 
with the same energy. We have obtained the transition 
line T C (H) for both Q1D and Q2D cases. As H becomes 
large, T C {H) increases oscillationally in both cases. 
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Fig. 3. Effective coupling constant as a function of aG/(27r) in 
the case of t c /t a = 0.02 and T/t a = 0.001. 



The reason why g/go is small for 4/4 = 0.3 is as fol- 
lows. In Fig. 4, we plot the effective coupling constant 
as a function of 4/4- In the case of 4 = 0, there is 
a logarithmic divergence at about 4/4 ~ 0.3, which is 
van Hove singularity for the quarter-filled band. 
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Fig. 4. Effective coupling constant as a function of t^/ta at 
aG/(27r)=0.025 and T/ta = 0.001. 



Therefore, the effective coupling constant normalizes by 
that for 4 = is small for 4/4 ~ 0.3. 

In this paper we have neglected the Zeeman term for 
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